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Starting with the general description of a moving curve, we have recently pre- 
sented a unified formalism to show that three distinct space curve evolutions 
can be identified with a given integrable equation. Applying this to the non- 
linear Schrddinger equation (NLS), we find three sets of coupled equations 
for the evolution of the curvature and the torsion, one set for each moving 
curve. The first set is given by the usual Da Rios-Betchov equations. The 
velocity at each point of the curve that corresponds to this set is well known 
to be a local expression in the curve variables. In contrast, the velocities of 
the other two curves are shown to be nonlocal expressions. Each of the three 
curves is endowed with a corresponding infinite set of geometric constraints. 
These moving space curves are found by using their connection with the in- 
tegrable Landau- Lifshitz equation. The three evolving curves corresponding 
to the envelope soliton solution of the NLS are presented and their behaviors 
compared. 
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1 Introduction 



The study of possible links between intrinsic kinematics of space curves 
and integrable soliton bearing equations p| deserves attention because of a 
wide variety of applications of moving curves such as vortex filament motion 
in fluids 0, dynamics of continuum spin chains 0, scroll waves in chemical 
reactions superfiuids 0, interface dynamics [0, etc. The pioneering work 
by Hasimoto on the motion of a vortex filament in a fluid was the first to 
suggest such a link. In classical differential geometry, a space curve embedded 
in three dimensions is described using the following Frenet-Serret equations 
for the ortho normal triad of unit vectors made up of its tangent t, normal 
n and the binormal b: 

ts = Kn ; Us = —Kt + rb ; = — rn. (1-1) 

Here, s denotes the arclength. The parameters K and r represent the cur- 
vature and torsion of the space curve. For a moving curve, these will be 
functions of both s and time u. The subscript s denotes partial derivative. 

Using the so-called local induction approximation. Da Rios [|] showed 
that the curve velocity v(s, u) at each point s of a vortex filament (regarded 
as a moving space curve) is given by the local induction equation : 

v = r„= Kh (1.2) 

For non-stretching curves, using compatibility conditions Tus = fsu = = 
Kgh — Km, along with similar compatibility conditions for t and n leads to 

i 

K^ = -{Kt)s-K,t; T^ = [{KjK)-r%-KKs. (1.3) 

Subsequently, the above coupled equations for K and r were also derived 
independently by Betchov , and hence Eq. ( p..3| ) is referred to as Da Rios- 
Betchov (DB) equations. Connection between these curve equations and 
solitons was unravelled by Hasimoto 0, whose analysis essentially showed 
that by using the following 'Hasimoto function' 

ip = K exp[i / r ds], (1.4) 
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the two equations in ( |1.3| ) can in fact be combined to give the integrable, 
sohton-bearing nonhnear Schrodinger equation (NLS) 

iqu + g.s + ^ \q? q = o, (1.5) 



with q = tp. Motivated by this resuh, Lamb [|T0[ gave a general procedure 
which helps identify a certain space curve evolution with a given integrable 
equation. Examples such as the NLS, the sine-Gordon equation and the mod- 



ified KdV equation were considered. Recently, we presented ITlJ a unified 
analysis to show that in general, two more distinct curve evolutions can also 
be so identified. In this paper, after outlining this analysis, we specialize to 
the NLS, and focus on the intrinsic kinematics of the three moving curves 
associated with it. We obtain coupled evolution equations for the curvature 
and torsion of each of the two new curves. These are certain analogs of the 
DB equations ( |1.3D . Each of the three moving curves is shown to be endowed 
with an infinite set of geometric invariants. Their natural connection with 
the integrable Landau-Lifshitz equation is pointed out, and a procedure to 
use its solution to construct the position vectors generating the three curves 
associated with the NLS is given. As an example, the moving curves for 
an envelope soliton solution of the NLS are found, and their behaviors com- 
pared. Interestingly, the velocities (at every point) of the two new moving 
curves underlying the general NLS evolution turn out to be certain nonlocal 
functions of the curve variables, quite in contrast to the local expression (Eq. 
(|1.2| )) for the velocity of the usual moving curve that has thus far been asso- 
ciated with the NLS. Possible application of these results to vortex filament 
motion in fluids is suggested. 

2 Identification of three curve evolutions with 
a given integrable equation 

To describe a moving curve, we find it convenient to introduce |12| the fol- 
lowing time evolution equations for the Frenet triad (t, n, b): 

t,, = gn + hh ; n„ = -gt + rob ; b„ = -ht - Ton. (2.1) 
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The parameters g, h and tq are, at this stage, general parameters which 
determine the time evolution of the curve. They are functions of both s and u. 
The subscript u stands for partial derivative. On imposing the compatibility 
conditions 

^SU ^US 1 ^SU ^US 1 ^su by^, (2.2) 



and using Eqs. (|LT|) and ( p.l|) , we obtain 



Ku = igs-rh); T^ = {TQ), + Kh] h, = {KTQ-Tg). (2.3) 



Lamb's procedure will be referred to as "formulation (I)", to distinguish 
it from two other additional formulations that are possible |TlI]. We shall use 
the subscripts 1,2 and 3 on the various curve parameters appearing in the 
three formulations, for ready reference. 

Formulation I: Here, one combines the second and third of Eqs. ( p..l|) , 
which immediately suggests the definition of a complex vector N = (n + 
ih) exp[2 /r (is], along with the complex function il) given in Eq. ( |1.4|) . By 
writing Ns,ts,N„ and tu in terms of t and N, imposing the compatibility 
condition N<ju = '^usi and equating the coefficients of t and N in it, one 
obtains 

i^u + lis + {l/2)[l (7i^*-7i» ds\i, = Q. (2.4) 



where 

7i = —{g + ih) exp[i j t ds]. (2.5) 

Thus for an appropriate choice of 71 as a function of ip and its derivatives, 
a known integrable equation for ip can be obtained from ( |2.4| ). Then, by 
comparing a known solution of this equation with the Hasimoto function 
(|1.4| ), the curvature K and torsion r of the corresponding moving space 



curve can be found. Next, using the above mentioned specific choice of 71 
in Eq. (p.5| ), the curve evolution parameters g and h can be determined as 
some specific functions of K ,r and their derivatives. Armed with these. 
To can also be found from the third equation in (|2.3| ). Thus a set of curve 
parameters K, t, g, h and tq has been found explicitly. In other words, a 
certain moving curve described by these parameters has thus been identified 
with a given solution of the integrable equation for ip. 



4 



As already mentioned, we have recently shown |]Tl|] that in addition to 
the above, there are two other distinct curves that can be identified with a 
given integrable equation. The procedures leading to these curves will be 
called formulations (II) and (III), respectively. 

Formulation (II): Combining the first two equations in (|1 . 1|) appropriately, 
we see that a complex vector M = (n — it) exp[i J K ds] and a complex 
function 

$(s,m) =rexp[i j K ds], (2.6) 

appear. Proceeding in a fashion analogous to (I) above, and imposing Ms„ = 

M.US, we get 

^u + l2s + {im[j (72$* -72*^) ds]^ = Q. (2.7) 

where 

72 = -(to - ih) exp[i j K ds]. (2.8) 

Formulation (III): Here, we combine the first and third equations of ( p..l| ), 
leading to the appearance of a complex vector P = (t — ib), and a complex 
function x given by 

x{s,u) = {K + it). (2.9) 
On requiring P^^, = P^^^, we get 



Xn + 73s + (l/2)[y (73X*-73X) ds]x = ^. (2.10) 

where 

lz = -{g + iro). (2.11) 
Now, Eqs. ( |2.7| ) and ( |2.10|) have the same form as Lamb's equation ( p.4|) . 



The discussion in the paragraph following Eq. ( p.5|) makes it clear that for 
an appropriate choice of 72 as a function of $ and its derivatives, and of 73 
as a function of x ^-iid its derivatives, known integrable equations for $ and 
X can be obtained from these equations. 

It is important to note from Eqs. ( |1.4| ), ( |2.6| ) and ( p.9| ) that the complex 
functions ■?/', $ and x that satisfy the (same) integrable equation in the three 
formulations are different functions of K and r. Further, we see that the 
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complex quantities 71, 72 and 73 (Eqs. (p.5D, ( p.8| ) and (|2.11| )) that arise in 
the three formulations respectively, also involve different combinations of the 
curve evolution parameters g, h and tq. Thus it is clear that the three for- 
mulations yield different sets of curve parameters, and hence describe three 
distinct moving space curves that can be associated with a given solution 
of an integrable equation. (We remark that it should be possible to extend 
our analysis to include partially integrable equations as well.) In the next 
section, we apply these results to the NLS, to demonstrate this clearly. 



3 The NLS : Two analogs of Da Rios-Betchov 
equations 

From the three formulations discussed in the last section, it can be easily 
verified that the respective choices 

71 = -ii>s; 72 = -i^s] 73 = -iXs, (3.1) 

when used in Eqs. (^), (^) and (gJ^) lead to the NLS given in (jT^), with 
q identified with the three complex functions ip, $ and respectively. Since 
the NLS is a completely integrable equation 0, it possesses an infinite set of 
conserved quantities J^, k = 1,2....cxd, which are in involution pairwise. The 
first three of these invariants are given by [|l^ 



A = y" \q\lds]l2= {l/2i) J {qsq* - qsQ*) ds] h = j - 



2] ds-.... 
(3.2) 



Let us consider the implications of this in the three formulations: 
(I) Setting q = if) in (O) and equating real and imaginary parts leads to DB 
equations (|1.3|). Therefore these also possess an infinite number of conserved 
quantities, obtained by setting q = if) m Eq. (^). The invariants now 
appear as geometric constraints [ll4| : 



h = j ds- l2 = I K^T ds; h = j [Kl + K^T^~ ^K^] ds;.... (3.3) 

Further, an inspection of Eqs. (|1.3|) shows that the total torsion Iq = J t ds 
is also conserved. This additional geometric constraint has no counterpart 
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among the NLS invariants given in ( p.2|) . 

(II) Here, we set g = $ in ( |1.5| ) to yield the coupled equations 

Ku = [{tsJt) - K% - TTs; Tu = -{Kt)s - TsK. (3.4) 



This is the first analog of the DB equations ( |1.3| ). As is obvious from com- 
paring ip in ( |1.4| ) and $ in Eq.( |2.6| ), this analog may be found by simply 
interchanging K and r in ( p..3| ). Thus the associated infinite number of geo- 
metric constraints can also be found using this interchange in Eq. (|3T 



Ii = J ds; l2 = J K ds; h = j [t^ + t^K^- ^r^] c/s; .... (3.5) 

Here, the total curvature Iq = J K ds is also conserved. This has no coun- 
terpart among the conserved densities ( |3.2| ) of the NLS equation. 
(Ill) Finally, setting g = x in ( |1.5| ) yields the following second analog of the 
DB equations: 

K^ = -T,,-^{K^ + T^)r; r.^ = K,^ + ^{K^ + t')K. (3.6) 



Next, setting g = x in Eq. (|3.2|) , we get the third set of infinite geometric 
constraints: 

h=J {K' +T^) ds; h = j {K,r -Kr,) ds; h = j [KI+tI~{K^+t'')] ds; .. 

(3.7) 

We end this section with the remark that the total length of the curve, 
L = J ds is also conserved in all three formulations, since the curves are 
non-stretching. 

4 Construction of the three curves using 
the Landau-Lifshitz equation 



A general solution of the NLS (Eq. ( |1.5| )) is of the form q = pexp[z^^]. 
Equating this with the three complex functions defined in Eqs. (|1.4| ), ( p.6|) 
and ( p.9|) will yield the curvature and the torsion of the three space curves to 
be (I) Ki = p, Ti = 9s, (II) K,2 = ds, T2 = p and (III) ^3 = pcos6, T3 = 
psin6. These are clearly, three distinct space curves, each with a known 
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curvature and torsion. However, solving the Frenet-Serret equations ( |1 . 1| ) 
explicitly for t by using these, to subsequently construct the position vector 
r(s, u) = ft ds of the (evolving) space curve, is a nontrivial task in practice. 
In the present context, we shall show that there is a certain link between 
the three curve evolutions and the integrable Landau Lifshitz (LL) equation 



T7| , which provides us with an alternative procedure to construct the three 
moving curves. 

We proceed by first equating the expressions for 71, 72 and 73 given in 
Eq. Q with those given in Eqs. (U), Q and (|2ll|) . This yields the 
following curve evolution parameters g, h and tq in the three cases: 

(I) gi = -KiTi; hi = Ku; Tqi = {nus/f^i) - r^, (4.1) 

(II) g2 = (rass/ra) - kI; h2 = -T2s\ = -^,2^2, (4.2) 

(III) (?3 = -r3.; /i3 = (l/2)(/€^ + r|); r^^ = k.^- (4.3) 
Substituting these expressions for each of the formulations appropriately in 



Eqs. (|2.1|), and using Eq. (|1.1[), a short calculation |Tl| shows that the 



completely integrable LL equation 



S X S,,; = 1 (4.4) 



is obtained in every case, i.e., Eq. ([4.4|) is satisfied by the tangent ti of 



the moving space curve in the first formulation, by the hinormal b2 in the 



second, and the normal in the third ||18[. Of the above, the first may be 
regarded as the converse of Lakshmanan's mapping where, starting with 
the LL equation, and identifying S with the tangent to a moving curve, one 
obtained the DB equations, and from them, the NLS for %p. The converses 
of the other two will yield the two analogs of the DB equations obtained in 
Section 3, and clearly lead to new geometries connected with the NLS. 
The LL equation ( ^4.4| ) has been shown to be completely integrable |T9 



and gauge equivalent EO] to the NLS. Its exact solutions can be found |19, 21 



We now outline how ri, V2 and r3, the position vectors generating the 
three moving curves underlying the NLS, can be found in terms of an exact 
solution S of Eq. (|4.4|) . Let (tj,nj,bj), i = 1,2 and 3, respectively, denote 
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the Frenet triads of the three moving curves satisfying Eqs. ( |1 . 1| ) , with the 
corresponding parameters and Xj. 

(I) Here, ti = ri^ is the tangent to a certain moving curve created by a posi- 
tion vector ri(s, u). Here, we set ti = S, a known solution of Eq. ( [4.4|) .Thus 

Ki = \tu\ = \Ss\; n = -2 = — (4.5) 

Thus the underlying moving curve ri{s,u) is given by 

ri{s,u) = j ti ds = j S(s,m) ds (4.6) 
The above expression for ri coincides with the surface that one would obtain 



using Sym's Jig soliton-surface method. 

(II) Here, b2 is the binormal of some moving curve r2(s, u). We set b2 = S. 
We have, t2 = r2s. The curvature K2 = ri and torsion T2 = ki (See Eq.(^]5|)). 
The position vector r2(s, u) generating the moving curve can be shown to be 



22 



r2{s,u) = [ t2 ds = / X — ds (4.7) 

J J Ki 

(III) Here, the normal of some other moving curve r3{s,u). Thus we set 
n3 = S. Further, = r^^. Next, from the Frenet-Serret equations ( |1 . 1[ ) for 
this triad, 

{kI + T^)ts = Ts{n3 X nss) - ksHss (4.8) 



A short calculation shows ||2^ that = Kicosrji and = ki sinr^i, where 
?7i = [/ti ds + Ciiu)]. Here, Ci{u) is a function of time m, which can 
be determined in terms of ki and ri, after a short calculation, by using the 
equations for ^^.u and t^u obtainable from Eqs. ([2.31). Then, substituting the 
values for and into Eq. ([4. 81), and setting na = S, the position vector 



r3(s, u) creating the third moving space curve can be found to be 
r3(., u)=Jtsds = J [(SxS^)^in,,-S. cos,,]) 
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5 Moving space curves associated with the 
NLS sohton. 



As seen in the last section, the expressions for the three moving curves asso- 
ciated with the NLS are given in (|4.6| ), (|4.7| ) and ([4 .91) respectively, in terms 
of a known solution S of the LL equation (|4.6|) . By defining three orthogonal 
unit vectors ei = |l,0,o|, §2 = |o, cos r], sin r^j, §3 = |o, — sin 77, cos r/j, a 
soliton solution |^ of ( [4 .61) can be written in the form 

S{s,u) = (1 — /iz/sech^(z/^))ei+/ii/sech(z/^) tanh(z/^)e2 — /iAsech(i/^)e3, (5.1) 

where, ^ = (s — 2Am), t] = (As + (z/^ — A^)m) and fi = 2i//(z/^ + A^). Here, u and 
A denote arbitrary constants. As is clear from our discussion in Section 4, 
the three curves that result from this solution are associated with the soliton 
solution q = 2z/sech(z/,^) exp {irj) of the NLS (Eq. (|1.5| )). They are found as 
follows: 

(I) Substituting Eq. (|5.1| ) in Eq. (^I5|), we get ki = 2z/sech(z/,^) and ri = A. 
Next, substituting Eq. (|0|) in Eq. (|4^ ) yields 



Ti = (^s — /itanh(z/,^), — yusech(z/^) cosrj, — /isech(z/,^) sin (5.2) 



This is seen to agree with the result obtained in |T5[ using Sym's ||T6[ pro- 
cedure. In Fig. (1), we have presented a stroboscopic plot of the moving 
curve ( |5.2| ) at different instants of time. This describes the propagation of 
the well known Hasimoto "loop" soliton. In between the times at which we 
have plotted the curve, the loop changes its size and also rotates about the 
asymptotic direction which the curve takes on at s ^ ± 00. In the figure, we 
have not plotted these intermediate times for the sake of clarity. The loop 
regains its form at subsequent times, as we see from the figure. 
(II) Here, we get ^2 = A and T2 = 2z/sech(z/^). Substituting Eq. i ^A\j in 
Eq. ( [4.7|) , the moving curve is found to be 

r2 = / [yuAsech(z/^)ei - /lA tanh(z/^)e2 + --15 ^§3] ds (5.3) 

In this case, since the curvature is constant and the torsion vanishes as s — >■ 
±00, for all finite times, the curve is bounded by two planar circles at both 



10 



ends. The maximum non-planarity of the curve occurs at s = 2 Am. A 
stroboscopic plot of Eq. ( p.3| ) given in Fig. (2). This curve is clearly seen to 
rotate and propagate as time progresses. 

(Ill) Here, ^3 = 2z/sech(i^^) cos r/ and T3 = 2z/sech(z/^) sin//. When Eq. 
(|5.1| ) is substituted in Eq. (|4.9| ), with rji = r], a long but straightforward 
calculation yields 



r3 = y [/isech(z/,^)(A sin 17 — z/tanh(z/^) cos 17)61 



— (/iA tanh(z/^) sin 77 + (1 — /^i^tanh (ly^)) cos 77)62 
(A^ — 1^"^) 

+ (tt^^ ^ sin 1] - fiX tanh(z/^) cos 77)63] ds (5.4) 

(A^ + u'^) 



A stroboscopic plot of the moving curve (|5.4|) at different instants of time is 
given in Fig. (3). This depicts the propagation of a "loop" soliton distinct 
from the Hasimoto soliton in the sense that it also oscillates with time, as is 
clear from the figure. Furthermore, at intermediate times, the loop rotates 
about the asymptotic direction, decreasing its loop size, and almost straight- 
ens out periodically, after which it starts looping again. These intermediate 
plots have not been presented in the figure. 

6 Three curve velocities associated with NLS 
evolution 

Returning to our general results of Section 4, the three position vectors rj, 
i = 1,2 and 3, which are associated with the NLS evolution, can be found 
from Eqs. ( [4.6|) , ( [4.7|) and (|4.9| ), respectively, using an exact solutions of the 



LL equation (|4.4|) . The corresponding curve velocities Vj(s,M) = rj„ at each 
point s, can therefore be found from these equations by direct differentiation 
with respect to time, u. However, to compare and contrast the intrinsic ge- 
ometries of the three space curves, it is instructive to express these velocities 
in terms of the vectors of the corresponding Frenet triads, as well as the 
curve parameters, as follows. Since the curves are non-stretching, we have 
Vjs = ^ius = Yisu = iiu- On the other hand, from the first equation of Eqs. 
(|2.1|), we have, tj„ = Qiiii + /ijbj. Here, the quantities Qi and hi for the three 
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curves of the NLS are given in Eqs. ( |4.1| ), ( [4 .21 ) and (^^), respectively. Using 
these in Vj(s,'u) = / tj„ ds, we get 



Vi = y ( — KiTiIli + Kisbi) (is = Ki hi. (6.1) 
V2 = y" [(r2ss/r2- /t2 ) n2 - K2sb2] c?s. (6.2) 

) V3= 1 [-r3,n3+ ^(/s:2 + r2)b3] rfs. (6.3) 



(I) 
(11) 
(III 

As expected, vi coincides with the vortex filament velocity Eq. (|1.2| ) derived 
by Da Rios in fluid mechanics. It is a local expression in the curve variables: 
The velocity vi at a point s depends on the curvature and binormal at 
that point only. In contrast, V2 and V3 are seen to be nonlocal in the curve 



variables: If we partially integrate the right hand sides of Eqs. (6^) and ( |6.3|) , 
and use the Frenet-Serret equations (Eqs. ( |1.1| )) repeatedly in the resulting 
expressions, both these velocities take on the form Vj = tj+ Bi nj+ Cj bj, 
i = 2,3, where the components Ai,Bi and Cj can be written in terms of 
integrals of certain functions involving the curvature, the torsion , their higher 
derivatives, and their various products. It is indeed interesting that in spite 
of such a complicated behaviour of these velocities, their corresponding curve 
evolutions are also endowed with an infinite number of constants of motion. 
This essentially stems from their connection with the NLS. More specifically, 
this is because, as we have shown in Section 3, the binormal b2 of the second 
curve and the normal of the third, satisfy the integrable LL equation. 

We conclude with the following remarks: In a fluid, as is well known, the 
induced velocity v at a point is determined as a volume integral involving the 
vorticity f2, by using the Biot-Savart formula. First, it is to be noted that in 
this formula, if one expresses fl in terms of the vectors of the Frenet triad of 
the filament, then, in any realistic model of a fluid, v is nonlocal in the curve 
variables, and becomes local only under certain approximations. Secondly, 
it is worth noting that in an interesting experiment with a fluid in a rotating 
tank, Hopfinger and Browand |^ had observed compact distortions which 
twist and propagate along a vortex core like a soliton. These two features 
suggest that our new results which link the velocities V2 and V3 with the 
NLS and its soliton, may be of some relevance in actual fluids. In view of 
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this, appropriate theoretical modehng of the vorticity, as well as possible 
experimental studies of the detailed geometric structure of moving vortex 
filaments, to look for any similarities with Figures (1) to (3) (which depict 
propagation of compact distortions) would indeed be of interest. 
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Figure 1: A stroboscopic plot of the evolving space curve ri(s, u) (Eq. ( |5.2| )) 
for V = 1 and A = 0.1. 




Figure 2: A stroboscopic plot of the evolving space curve r2(s,M) (Eq. (|5.3| )) 
for V = 0.3 and A = 0.1. 




Figure 3: A stroboscopic plot of the evolving space curve r3{s,u) (Eq. (|5.4| )) 
for V = 1 and A = 0.3. 
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